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Abstract 

We discuss forward-backward correlations in the multiplicity of 
produced particles in heavy ion collisions. We find the Color Glass 
Condensate generates distinctive predictions for the long range com- 
ponent of this correlation. In particular, we predict the growth of the 
long range correlation with the centrality of the collision. We argue 
that the correlation for baryons is less strong than that for mesons. 



1 Introduction 

The Color Glass Condensate provides a phenomenologically successful de- 
scription of high energy hadronic processes [1]-[11]. The basis of the descrip- 
tion for the infinite momentum hadronic wavefunction is to decompose the 
wavefunction into pieces associated with the fast moving constituents, and 
a piece associated with the slower moving components [9]. The fast moving 
components act as sources for the slow moving components. The slow moving 
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components are treated as classical fields. Because the density of the gluons 
in the hadronic wavefunction is large, the coupling constant is weak. The 
field strengths are large, nevertheless, A ~ 1/g, and although one can use 
weak coupling techniques, the problem is essentially non-perturbative due to 
the strong field. 

The physical picture of the soft component of the hadronic wavefunction 
can be easily understood in terms of the phase space distribution 



dyd'^PTd'^rT 

For pt < Qsat, where Qsat is the saturation momentum, rj ~ l/a^. This 
means that the quantum mechanical gluon states are multiply occupied, and 
saturate at l/ag. At this phase space density, repulsive interactions become 
large. As one adds more gluons to the system, these components of the 
hadronic wavefunction remain fixed, since the repulsive interaction energy 
makes this unfavorable relative to adding in a gluon in a less occupied state 
at higher transverse momentum. Above the saturation momentum, the phase 
space density becomes small. When one adds more gluons to the system, they 
typically have momenta above and near the saturation momentum. This in- 
creases the saturation momentum. The renormalization group equations for 
the Color Glass Condensate predict that the saturation momentum never 
saturates, that is, it grows with decreasing x of the gluons in the wavefunc- 
tion [12]- [15]. One is perpetually adding in more gluons at an ever increasing 
saturation momentum. 

When applied to hadron-hadron collisions, one imagines the collision of 
two sheets of colored glass [16]-[23]. Initially, the CGC fields arc transversely 
polarized color electric and magnetic fields. At a very short time after the 
sheets collide, t ~ e"''/"", the fields become longitudinal electric and magnetic 
fields in the central region between the colliding nuclei, and the transverse 
field in this region vanishes. The original transversely polarized fields are 
associated with the fast moving components of the nuclei and remain intact. 
This situation is shown in Fig. 1. The matter produced immediately after the 
collision is called the Glasma. It has properties distinct from the Color Glass 
Condensate and the Quark Gluon Plasma. It exist at times intermediate 
between the two, and has some properties in common with both the CGC 
and the QGP, hence the name. 
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Figure 1: (a) Two sheets of colored glass approaching one another, (b) After 
the collision, an equal and opposite charge density of color electric and color 
magnetic charges are set up on each nucleus. This induces longitudinal color 
electric and magnetic fields in the region between the nuclei. 

The origin of the longitudinal fields arises because as the sheets of colored 
glass pass through one another, the fast moving components have added to 
them a distribution of color electric and color magnetic charge. The charge 
density in the transverse plane of one nucleus is the negative of that in the 
other nucleus. 

This picture of the fields has a restricted range of validity. For the effective 
action which describes these fields to be valid, we require that the field exist 
over a region of rapidity 



In a larger range of rapidity, quantum fluctuations become important, and 
the description in terms of classical fields breaks down. 

It has also been argued that such long range correlations exist in the dual 
parton string picture of hadron collisions [24, 25]. Some time ago, it was 
also argued that the Color Glass Condensate has long range correlations in 
rapidity [26] . The Glasma provides a theoretical framework following directly 
from QCD which includes many features of the dual parton model approach, 
in particular the longitudinal rapidity structure. 

In this paper, we explore these long range correlations. We use the for- 
malism of forward-backward correlations. We argue that there is a long range 



Ay < l/a, 



(2) 
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component of fluctuations for the gluons, 



/dNf}dNlj\ 1 dm^ 
\ dyi dy2 / dy ' 

which is approximately rapidity independent. Here F (B) are regions in 
the positive (negative) rapidity region, and y is some rapidity intermediate 
between yi and y2, oftentimes taken to be the center of mass rapidity. For 
quarks, 

/d_Nldn\d_m_ 
\ dyi dy2 / dy 

The essential difference between the two is that for the gluons the is a factor 
of l/ttjj. This can be understood since the gluons are acting coherently. The 
decay process involves the coherent emission of l/cts gluons. In the language 
of clusters, the typical cluster size of produced by a fireball is 1/as particles. 
The quarks are fermions, and their typical cluster size is of order one, since 
there is no coherence for fermions. 

To compute the forward-backward correlation function, we also need to 
know the short range component. We define^ 

_ {NpNb) - {Nf) {Nb) 
- {m) - {Nf ' 

with N the multiplicity in the forward or backward interval. We argue below 
that there is a short range component which falls exponentially in rapidity. 
It is of strength of order relative to that of the long range component. 

We will argue below that at short range, the correlation function is pol- 
luted by final state effects, and is not directly computed from the initial 
state distribution of gluons. If we look at the long range piece of (Tfb, say 
at rapidity separation of greater that one unit of rapidity, final state effects 
should not be so important. 

We argue below that this long range component has the form 

ctfb = ^ 2 ■ (6) 
1 + cat 



^This quantity coincides with the h usually defined through {nF){nB) = a + briB [27]. 
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We do not perform a detailed computation of the constant c. It involves 
knowing the strength of the soft correlated emission. There is however a 
qualitative effect which we can understand about c: If the contribution of 
baryons to the correlation function is large, then c increases. This is because 
this contribution brings in no factor of a^. The contributions of baryons may 
be important of central collisions at large Pt- 

Notice that as a function of increasing centrality, decreases. The 
forward-backward correlation function should therefore increase for more cen- 
tral collisions. 



2 Computing the Forward-Backward Corre- 
lation 

Consider the collision of two nuclei in the center of mass frame. We treat 
the high momentum degrees as sources. There are sources for each nucleus. 
Let us consider some region of transverse extent a. The value of a will be of 
the order of the inverse transverse momentum scale which we probe. For the 
total multiplicity fluctuations, this size is of the order a ~ 1/Qsat- 

Now in the transverse region of size a, due to fluctuations, one of the nuclei 
will have a larger charge than the other. The situation becomes like that of 
pA collisions where there is an asymmetry between target and projectile [28]- 
[29]. The total multiplicity is given by the number of gluons in the object 
with the smaller charge, evaluated at the saturation momentum scale of the 
larger charged object. This is because the larger charge object is effectively a 
black disk for the smaller charged object, and as such, during the scattering 
the smaller charged object materializes all of its gluons. Since the difference 
in charge between these objects should be small compared to the charge (so 
we can neglect the difference between the Qsat oi projectile and target), for 
each region of size a, the total multiplicity is therefore of order 

a g,„, . 7 

ay as 

Summing over all areas gives the familiar Kharzeev-Nardi formula for the 
total multiplicity [30] , 
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Figure 2: The leading order diagram which induce long range correlation in 
rapidity. The source of one nucleus is given by the x and that of the other 
by the o. The produced gluon is denoted by the dotted line. 

The fluctuations are given by the square of the diagram for the multiplic- 
ity. The leading order contribution from the classical fields comes from the 
diagram of Fig. 2. (This diagram at first sight appears to be uncorrelated, 
but becomes correlated after averaging over the sources associated with the 
nuclei. This averaging "ties together" the sources in the upper and lower part 
of the diagram.) There is a factor of l/a^ which is due to the fact that this 
is the square of two classical processes. The correlation will only be non-zero 
if the same transverse area is probed. So the factor of a^, upon summing 
over areas becomes of order t^R^o^. There is also factor of Q^sat- ^^"^ total 
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multiplicity squared, a ~ 1/Qsat, this gives 




dN 



) 



2 



> 



1 dN 

as dy ' 



dy 



sat 



(9) 



The fluctuation at different rapidities is the same since the fleld which pro- 
duces these particles is rapidity invariant: 



Note the factor of l/ag in this result. This is because the l/cts gluons 
associated with the classical field behave coherently. Its origin is similar 
to that of the dependence on cluster size for fireball models of correlations. 
One needs to know the number of particles in each fireball decay. This 
number scales the square root in A^(oc l/cts) fluctuation since one assumes 
each flreball is produced statistically but an extra factor of N comes because 
each flreball decays into N particles. 

There is of course a correlated piece, which comes form the diagram of 
Fig. 3. Here one of the gluons is associated with the forward rapidity and 
the other with the backward rapidity. This diagram has two factors of ctg, 
and should give a contribution to the total multiplicity fluctuations of order 



Here k is a constant of order one. Note that dNcor/ dy is of order relative to 
that for the classical contribution. There is no interference with the leading 
order diagram of Fig. 2, since the average of an odd sources from the same 
nucleus must vanish. It should also be noted that outside of the context of 
perturbative QCD, there should be diagrams like Fig. 3 associated with the 
scattering of the hadronic degrees of freedom of the nuclei. 

To understand the different factors of ag between the diagrams of Fig 2 
and Fig 3, we can think of the source as strong, and therefore of order 1/g. 
Therefore each of the diagrams squared in Fig 2 is of order l/g'^, making an 
overall contribution of order l/ctg. In Fig. 3, the diagram is overall zeroth 
order in powers of g. 




(10) 




(11) 
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Figure 3: Here two gluon from the different sources scatter and produce two 
gluons in the forward and backward hemisphere. 



For the long range correlation function, we see that the correlated piece 
is overall a factor of smaller than the classical contribution. We therefore 
conclude that at large rapidities 



1 



1 



cat 



(12) 



where is to be evaluated at the saturation momentum associated the cen- 
trality of the collision and we have considered two intervals F, B symmetric 
around midrapidity. It clearly becomes more correlated as the centrality in- 
creases. It should also be noted that c is an increasing function of the rapidity 
difference Ay =| yi — y2 |- To see it, we define agC— ^^^/^ to get 



-kAj/ 



c — c 



1 + aice-^^y ' 



(13) 
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which clearly increases with increasing A|/. 

We can repeat all the steps above for baryons. The only difference is that 
the baryon multiplicity is of order a relative to that of the gluons, and that in 
the fluctuation formula, there is no additional factor of l/a^ associated with 
coherent emission. Thus the baryon correlation function should not grow 
with increasing centrality unlike of mesons which presumably ultimately arise 
from gluons. 

In the most central coUisions at RHIC, at high pt, there is considerable 
contribution from baryons, and so measurements of correlations at high px 
may not have such a large correlation for central collisions as would that for 

mesons. 

In general, the correlations as a function of pt may be interesting. The 
low Pt gluons are in a maximally occupied state. The repulsive interaction of 
the gluon in these states will not allow the addition of more gluons. Changing 
the charges on the nuclei therefore does not affect the multiplicity of low p^ 
gluons, so there should not be large correlation. The dominant contribution 
to the total multiplicity fluctuations comes from gluons with pt ~ Qsat- For 
Pt > Qsat, there should be considerable correlation. Of course the transverse 
momenta of gluons is strongly affected by final state interactions, and it is 
unclear how this effect could persist. 

As a final comment, let us recall the result for forward-backward corre- 
lations in a different framework: string models which allow for some kind of 
collectivity [31] -[33] e.g. a non-thermal phase transition like percolation. The 
generic finding in these models is a reduction of the long range correlations 
with increasing collectivity. While this reduction increases with increasing 
centrality, it does not imply that the correlation itself reduces with increasing 
centrality for a given modeling of collectivity. Indeed, if Ng is the number 
of independent particle sources and np the number of particles in the for- 
ward hemisphere from each source, the same assumptions leading to Eq. (12) 
result here in 

A is [34] a decreasing function with increasing centrality, so a pb increases, 
as in the CGC, with increasing centrality. 
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3 Summary and Conclusions 



The strongest objection to this analysis is that we have used quarks and 
gluons as degrees of freedom. In reaUty, we measure mesons and baryons. 
It is standard lore that the meson multiplicity distribution directly reflects 
that of the gluons. Baryons, particularly baryon-antibaryon pairs may be 
subject to more drastic final state interactions. Baryons can be enhanced at 
intermediate px in central collisions. 

For long range correlations, we have not imagined a final state process 
which can affect the correlation in the total multiplicity. At short range, 
diffusion can have a strong effect. Imagine that after particles have been 
produced, some of the particles from the forward region diffuse to the back- 
ward region. This decreases the number of forward particles and increase 
those in the backward direction. This can induce a negative contribution to 
the correlation function. Certainly for rapidity differences of less than one 
unit, processes such as this will be important. That said, understanding such 
diffusion and comparison with correlation function measurements may probe 
such effects, and may be interesting in their own right. 

Fluctuations in the correlation between centrality and impact parame- 
ter may give long range correlations. It is difficult to understand how the 
correlation could become stronger with more central collisions. Such corre- 
lations can be estimated by a variety of Monte Carlo event generators, and 
the strength of any resulting signal evaluated. 
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